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ABSTRACT 

We present high-resolution HIFI spectroscopy of the nucleus of the archetypical starburst galaxy M82. Six '^CO lines, 2 '^CO lines and 4 
fine-structure lines have been detected. Besides showing the effects of the overall velocity structure of the nuclear region, the line profiles also 
indicate the presence of multiple components with different optical depths, temperatures, and densities in the observing beam. The data have been 
interpreted using a grid of PDR models. It is found that the majority of the molecular gas is in low density (n = 10' cm"') clouds, with column 



densities of A'h = 10" cm 



X 
J3 



and a relatively low UV radiation field (Go= 

12.75 X 



10"). The remaining gas is predominantly found in clouds with higher 

10' ' 



dominated by a small (1% relative surface filling) component, with an even higher density (n= 10 cm ) and UV field (Go= 
show the strength of multi-component modelling for interpretating the integrated properties of galaxies. 
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- cm ^). The highest J CO lines are 



10 ). These results 



1. Introduction 

M 82 is one of the best studied starburst galaxies in the local uni- 
verse. Its short distance (3.9 Mpc, Sakai & Madore 1999) makes 
it a superb candidate for detailed studies of the physical pro- 
cesses related to star formation and their effects on the galaxy. 
The emission from the interstellar medium (ISM) is an excellent 
tool for diagnosing the physical and chemical properties of the 
star-forming environment. Over the past decades, M 82 has been 
studied in many atomic and molecular species. These observa- 
tions show a complex environment where multiple components 
with different exci tation, temperatu r es, densities, and filling fac- 
tors coexists (e.g. . Wild et al. 1992; Lord et al. 1996; Mao et al. 



20001; IWeiB etal] [2001: Ward et al.. 2003: .Spaans & Meiierinfc 



20071 iFuente et alj|2008h . 



In this paper we present observations of the nucleus of 
M82 using the Heterodyne Instrument for the Far Infrared 
(HIFI, De Graauw et al., 2010) on board of the ESA Herschel 
Space ObservatorjQ as part of the HEXGAL key programme (PI 
R. Glisten). Due to the large spectral coverage available with 
Herschel, we can observe a large number of lines, enabling com- 
prehensive study of the excitation of the different ISM compo- 
nents. At the same time, the high spectral resolution provided 
by HIFI makes it possible to separate different velocity compo- 
nents. In this paper, we combine these observations with detailed 
modelling to derive the physical conditions and excitation mech- 
anisms of the nuclear ISM of M 82. 



* e-mail: loenen@strw.leidenuniv.nl 

' Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



2. Observations, reduction, and results 

The observations were performed in two blocks during April 17- 
19 and May 03, 2010. Selected CO transitions and their '^CO 
isotopomers, as well as fine-structure lines of [CI], [CII], and 
[Nil], were observed (see Table [T] and Fig. [1]) in the nuclear 
region of M82 (RA 09h55m52.22s, Dec 69d40m46.9s J2000). 
The observations were carried out in dual-beam switch mode 
with standard 3' beam-chopping in the "fast" mode with chop 
rates between 0.2 and 2 Hz to improve stability. Relevant details 
of the observations are summarized in Table [T] 



2.1. Data reduction 

Initial data processing was done using HIP^ version /.9. A 
modified version of the level 2 pipeline was used that does not 
time average the data, in order to inspect the individual sub- 
scans in observations that contain more than one subscan. After 
inspection, the subscans were averaged and the subbands were 
stitched together. 

Further analysis of the data was done using the CLASS pack- 
age0. The data were binned until a sufficiently high signal to 
noise ratio was achieved. Linear baselines were subtracted, ex- 
cluding parts of the baseline that contained lines or seemed af- 
fected by instabilities. After this subtraction, the '■'CO (5 -4) and 
'^CO (6 - 5) spectra showed residual baseline structure. Since 
these structures are not seen in the '^CO (5 - 4) and '^CO (6 - 5) 
spectra, they were assumed to be baseline instabilities and were 
removed by fitting a 4* order polynomial. 



Herschel Interactive Processing Environment 
' Continuum and Line Analysis Single-dish Software: 
|http://www.iram.fr/IRAMFR/GILDAS. 
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Fig. 1. Observed spectra. Both polarizations were averaged in these spectra. Velocities are heliocentric. 



2.2. Spectra 

Figure [T] shows the resulting spectra, obtained by averaging both 
polarizations. Six CO lines, 2 '^CO lines, and all 4 targeted fine- 
structure lines are detected. Three lines are clear non-detections: 
CO (13 - 12), i^CO (9 - 8), and '^CO (10 - 9). The '^CO (7 - 6) 
line can be seen in the spectrum, but the integrated flux can not 
be securely determined due to severe baseline instabilities. The 
'^CO (8 - 7) observations also suffer from baseline problems be- 
tween ~ 400 and 700 km s"', making it impossible to determine 
a good baseline and establish the presence of a line. Because of 
the high frequency of the lines, the [C II] and [N II] spectra cover 
only a narrow range, and therefore a baseline can only be fit to 
the outer channels. 



suits in values for Ks between 1 and 0.19 (see Table [T]). Since 
the calibration of HIFI is still preliminary, theoretical predictions 
(Kramer'2006) are used to convert the antenna temperatures into 
flux densities (see Table[TJ. After applying the beam corrections 
and the preliminary calibration, the final CO fluxes are very sim- 
ilar (within 10%) to the fluxes measured with SPIRE, which 
are c orrected for source-b eam coupling using a 250 jum SPIRE 
map (Panuzz o et al.l 1201 0). The final corrected fluxes are listed 
in Table [1] together with the Icr uncertainties derived from the 
fits. Given the uncertainties in the pointing and calibration, we 
apply an additional error to the fluxes, when comparing them 
to the models. This error will increase with frequency; however, 
since the exact frequency dependence is unknown, an average 
value of 30% is used. 



2.3. Line fluxes 

The CO emission in M 82 has two main velocity components 
(the southwest (SW, Vhei ~160kms"') and northeast (NE, 
Vhei -300 kms^') lobes; see e.g.. Fig. |2] and Wild et al. 1992), 
which can also be seen in our data. Therefore, two Gaussian 
profiles are fitted to each spectrum. The size of the beam varies 
for the different frequencies, ranging from 12" for the [C II] ob- 
servations to 44" at the lowest frequency ([CI]609/im). To be 
able to compare the observations to each other and to the model 
results, all observations are scaled to the largest beam. The scal- 
ing factor (ks) is calculated by convolving a 450//m SCUBA 
map with a resolution of 7" (taken from the SCUBA archive) 
with Gaussian profiles with the size of the different beams, to 
estimate the flux contained in those beams. The value of a-, is 
calculated by taking the ratio of the integrated 450 fim flux for 
each beam with the flux contained in the largest beam. This re- 



3. Analysis and discussion 

The spectra in Fig.[T]show well-resolved line profiles, which pro- 
vide important boundary conditions for further analysis. Most 
lines show the two main velocity components: the blue-shifted 
SW and the red-shifted NE emission lobes. The [Nil] and [CII] 
lines show spectra that are dominated by the blue-shifted com- 
ponent. This can occur because the beam at these frequencies is 
so small (see beam size bars overlayed on Fig. |2]i that the NE 
lobe is not covered. Our observing position is shown in Fig. 3 
of WeiB et al. (this volume), where various beam sizes are also 
indicated, showing the extent to which the two lobes are covered 
at various frequencies. 

Inspection of the '^CO line profiles shows that at low J 
the SW lobe, peaking at 160 km s"', dominates the blue-shifted 
emission component. Going to higher J, an additional feature 
peaking at 100 km s"' becomes increasingly important and is 
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Fig. 2. Position-velocity diagram of the Nucleus of M 82, ob- 
served in '^CO(6 - 5). Vertical resolution and beam size are 
~16". The two red bars show the largest (44") and (12") small- 
est beams in our observations. 

brighter than the 160kms ' peak at / > 8. This feature dom- 
inates the '^CO lines already at 7 = 5 - 4. '^€0(3 - 2) and 
'^CO (6 - 5) position-velocity diagrams (Fig.|2] data taken from 
the JCMT archive) show that this 100 km s ' feature corresponds 
to a separate structure in the SW lobe. The change in contri- 
bution of this feature for different transitions and isotopomers 
indicates that it has physical properties different from the neigh- 
bouring 160 km s ' feature. This demonstrates that M82 con- 
tains star-forming regions with different physical properties, 
such as optical depth, temperature, and density, in the observ- 
ing beam. This information only becomes available with the 
velocity-resolved line profiles provided by HIFI. 

3.1. CO excitation 

The NE and SW lobes, unlike the lOOkms ' spectral feature, 
do not show any difference in behaviour. They merely reflect 
the kinematics of M 82 and will both have contributions from 
different physical components of the ISM. Because of this and 
the 100 km s"' component being too weak in most lines to be fit 
separately, we analyse the CO excitation using the integrated line 
strengths, but explicitly allowing for the presence of multiple 
excitation components. 

We combined our data with ground-based observations of 
CO ( 1 - 0) to CO (4 -3), CO (6 -5 ), CO (7 - 6), and '^CO ( 1 - 0) 
to I^CO (3 - 2) (IWard et alJl2003h and the SPIRE detections of 
CO (4 - 3) to CO (13 - 12) '^CO(5 - 4), '^CO(7 - 6), and 
i^CO(8 - 7) (iPanuzzo et aUlIoTol) . Imrdetan (HqoI present 
separate observations of the two lobes. Since both fit within our 
largest beam, we added the fluxes. The ground-based data have 
been scaled to match our ob servations , usi ng the CO (6 - 5) and 
CO (7 - 6) fines. Botfi the Ward et alJ (12003) and Panuz zo et all 
(12010) data have already been corrected for source-beam cou- 
pling factors. This resulted in the two CO excitation diagrams 
shown in Fig. [3] 

The CO and the fine-structure line ratios are compared to 
the resul ts of a large g rid of photon-do minated reg ion (PDR) 
models (iMeiierink & Spaans..2005; Meiierink et al.li200 7). The 
^^CO/^^CO abund ance ratio is fixed at a value of 40 (similar to 
the value found by 'Ward et al."2003V The comparison with the 
CO lines is shown in Fig. [3] The observations can be best re- 
produced with a combination of one low and two high-density 
components. The low-density component, which dominates in 
the lowest CO lines (J < 3), is found to have a density of 
n = 10^'^ cm"^ and a UV flux of Go- 10^. This component 
represents the extended molecular ISM. The mid-/ CO lines 
(4 < 7 < 7) are mostly produced by clouds that have a density 
of n = 10^ cm"-* and a radiation field of Go= 10^^^. At transi- 




upperj level 



Fig. 3. Excitation of '^CO and '^CO. Solid squares repre- 
sent HIFI observations presented in this work, and open dia- 
monds^ and circles data are taken from Ward et al. (2003) and 
Pan uzzo stall (l2010l) . respectively. Lines represent the result of 
PDR models. See legend and text for more information on the 
models. 

tions higher than 7 (and '^CO J > 5), the line emission comes 
from clouds with even higher densities (n - 10^ cm"^) and radi- 
ation fields (G()= 10^^^'). These two models represent the dense, 
star-forming molecular clouds. The densest component has pa- 
rameters similar to hot cores in the Milky Way. This component 
has excitation properties similar to the 100 km s"' feature de- 
tected in our spectra as can be seen by comparing the '^CO and 
'^CO (6 - 5) lines in the models. 

Since we have '^CO lines for both the low- (J < 3) and 
high-density (J > 4) components, it is possible to determine the 
column density for both. The low-density component is found 
to have a slightly higher column density (A^h = 10^'^ cm"^) 
than the high-density component (A^h = 10^''^ cm"^). This is 
reflected in the '^CO/'^CO line ratios, which increase from 14.5 
for the y = 1 - to 27.2 for the 7 = 8-7 transition. 

The relative scaling of the three models provides an estimate 
of the relative beam filling of the components and was found to 
be 70%:29%:1% for the low density, high density low Go, and 
high density high Go components, respectively. 

3.2. Fine-Structure lines 

To further constrain the modelling results, we also considered 
fine-structure line ratios. The ratio betwe en the two f C I] fines i s 
sensitive to the impinging UV flux (see IMeiierink et al.ll2b07h . 
Line ratios with [C II] and [N II] may provide additional insight 
into the ionization balance. Because of the different shapes of the 
red component of the [C II] and [N II] lines, we do not use the 
total integrated flux of those lines. Instead, the blue component 
is scaled using the average total/blue ratio. 

The model prediction for the [C I] 609 jum/[C I] 369 /vm line 
ratio (0.28) is close to the observed value (0.31). This is sur- 
prising since [C I] 609 fim is usually observed to be stronger 
than predicted by models or observed in the Milky Way (e.g., 
IWhite et alJll994l: llsrael & BaasI 120021) . That our prediction is 
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Table 1. Observation parameters and observed line fluxes, where fluxes include beam coiTection and calibration. 



Line 


Freq 


A 


Band 


tint 


Beam 




Umb 


iA-»Jy 


Blue 


Keu 


Total 




[Orlzj 






[sec] 


L J 






LJy ^ J 


[tu Jy Km s J 


[lU Jy Km s J 


[lU Jy Km s J 


CO y = 5 - 4 


576 


520 


lb 


250 


38 


0.87 


0.69 


474 


44.2+0.4 


34.3+0.4 


78.5+0.6 


CO y = 6 - 5 


691 


434 


2a 


248 


33 


0.77 


0.68 


481 


44.2±0.8 


32.5±0.8 


76.6+1.2 


CO y = 7 - 6 


807 


372 


3a 


248 


27 


0.61 


0.68 


481 


40.6±1.6 


30.5±1.7 


71.2+2.3 


CO y = 8 - 7 


922 


325 


3b 


250 


25 


0.55 


0.68 


481 


33.8±2.6 


26.4±2.7 


60.2+3.8 


CO y = 9 - 8 


1037 


289 


4a 


306 


23 


0.48 


0.67 


488 


31.4±1.4 


18.0±1.4 


49.3+2.0 


CO y = 10 - 9 


1152 


260 


5a 


334 


20 


0.40 


0.67 


488 


16.9±3.7 


19.5±4.3 


36.4+5.7 


'^"coy = 5-4 


551 


544 


la 


314 


44 


1.00 


0.69 


474 


2.0±0.1 


1.3±0.1 


3.4+0.2 


'^coy = 6-5 


661 


453 


2a 


306 


33 


0.77 


0.68 


481 


2.3±0.3 


0.8±0.2 


3.1+0.4 


[CIl^Pi -Spo 


492 


609 


la 


248 


44 


1.00 


0.69 


474 


13.3+0.6 


8.2±0.6 


21.5+0.8 


[CI] Pi 


809 


369 


3a 


306 


27 


0.61 


0.68 


481 


22.2±1.9 


20.3±2.1 


42.5+2.8 


[CII]2p3/2 -2 Pi/^ 


1901 


158 


7b 


306 


12 


0.19 


0.63 


519 


1051.4+57.6 


1262.7±57.0 


2314.1+81.1 


[NII]3Pi 


1461 


205 


6a 


507 


15 


0.27 


0.65 


503 


58.0±39.1 


66.7±39.1 


124.7+55.3 



close to the observed value stems from our multi-component 
modelling. The low-density component has a ratio that is higher 
than observed, and the high-density components have lower ra- 
tios, resulting in an average that is close to the observed value. 

The models under-predict the [C II]/[C I] ratios by about a 
factor of three higher. This excess in [C II] emission likely arises 
from the diffuse, ionized ISM, which is confirmed by the simi- 
larity of the [C II] and [N II] line shapes of the blue lobe. 

To address the uniqueness of our model, a model with only 
one high-density component was also attempted. The best fit is 
found using a lower density (n = 10''^'' cm"^) and a high radia- 
tion field (Go- lO''^). Also a much higher surface ratio is needed 
(ratio of low- to high-density components 20%:80%). Although 
this model is able to fit the '^CO excitation within the error bars, 
it underpredicts the strength of the '^CO lines by factors of a 
few. Also the predicted [C II]/[C I] ratios are an order of mag- 
nitude too high and, therefore, a two c omponent interpretation, 
such as done by iPanuzzo et al.l (l2010h . is insufficient for mod- 
elling the CO excitation in M82. This conclusion is reinforced 
by earlier analyses by e.g., based on large velocity gradien t anal- 
ysis of multiple CO lines (e.g.. Wild et al. 1992; Gust en et al.l 
119931: IWeiB et al.ll200ll IWard et al...2003^ .WeiB et al..,20Ql 

4. Summary and outlook 

The CO and fine-structure lines provide an excellent tool for de- 
termining the physical parameters of the nuclear ISM of M 82. 
The line shapes reveal distinct velocity components at 100, 160, 
and 300 km s"^ While the 160 and 300 km s"' features reflect 
the large-scale kinematics of M82, the 100 kms' feature shows 
an excitation different from the other components, which is only 
revealed by the velocity-resolved HIFI data. The presence of 
multiple ISM components with different physical conditions is 
also found in modelling of the CO lines, which shows that three 
density components are needed to explain the observed line 
fluxes. 

The majority of the gas is in low density (n = 10^'^ cm""*) 
clouds, with column densities of A^h - 10^''^ cm"^ and a rel- 
atively low UV radiation field {Gq= 10^). The remaining gas 
is predominantly found in clouds with higher densities (n - 
10^ cm"^) and radiation fields (Gq- 10^^^), but lower column 
densities (A'h = lO^^'^ cm"^). The highest J CO lines are domi- 
nated by a small (1% relative surface filling) component, with an 
even higher density (n = 10'' cm^^) and UV field (Go= lO^-^^). 

Using the unique spectral coverage and resolution offered 
by HIFI, we are one step closer to correctly dissecting the intri- 



cate and interacting density components and excitation regimes 
present in the core of M 82. Additional observations (part of the 
HEXGAL programme) of the lobes of M 82 and of other star- 
burst galaxies will increase our understanding of such systems 
even more. 
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